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ABSTRACT 

The analysis of relativistically broadened X-ray spectral features from the inner accretion disk 
provides a powerful tool for measuring the spin of supermassive black holes (SMBH) in active galactic 
nuclei (AGN). However, AGN spectra are often complex and careful analysis employing appropriate 
and self-consistent models are required if one is to obtain robust results. In this paper, we revisit the 
deep July-2009 Suzaku observation of the Seyfert galaxy NGC 3783 in order to study in a rigorous 
manner the robustness of the inferred black hole spin parameter. Using Monte Carlo Markov Chain 
(MCMC) techniques, we identify a (partial) modeling degeneracy between the iron abundance of the 
disk and the black hole spin parameter. We show that the data for NGC 3783 strongly require both 
supersolar iron abundance (Z-p c = 2 — 4Z Q ) and a rapidly spinning black hole (a > 0.88). We discuss 
various astrophysical considerations that can affect the measured abundance. We note that, while the 
abundance enhancement inferred in NGC 3783 is modest, the X-ray analysis of some other objects has 
found extreme iron abundances. We introduce the hypothesis that the radiative levitation of iron ions 
in the innermost regions of radiation-dominated AGN disks can enhance the photospheric abundance 
of iron. We show that radiative levitation is a plausible mechanism in the very inner regions of high 
accretion rate AGN disks. 

Subject headings: galaxies: individual(NGC 3783) - X-rays: galaxies - galaxies: nuclei - galaxies: 
Seyfert -black hole physics 



1. INTRODUCTION 

The X-ray spectra of active galactic nuclei (AGN) often 
display relativistically broadened spectral features from 
the innermost regions of the black hole accretion disk. 
Analysis of these features, most notably the broad iron 
line, is one of the most powerful tools that we currently 
possess for studying the region close to the supermassive 
black hole (SMBH) where the gravitational field is rel- 
ativistic (Fabian et al. 1995; Reynolds & Nowak 2003; 
Miller 2007). As such, essentially all robust measure- 
ments of SMBH spin are currently derived from an anal- 
ysis of these spectral features (Dabrowski et al. 1997; 
Brenneman & Reynolds 2006; Reynolds et al. 2012). 

However, the X-ray spectra of AGN are often com- 
plex and careful modeling is required if one is to obtain 
robust results. The broad iron line results from fluores- 
cence in the surface layers of the optically-thick accretion 
disk when it is irradiated by an external source of hard 
X-rays (the corona or the base of a jet). When using 
high signal-to-noise spectra to study black hole spin, it 
is insufficient to study the broad iron line in isolation; 
one must self-consistently model the full spectrum of X- 
rays that are reflected/reprocessed by the accretion disk 
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(George & Fabian 1991). Photoionization of the accre- 
tion disk surface must be modeled (Ballantyne, Fabian & 
Ross 2002; Ross & Fabian 2005; Ballantyne, McDuffie & 
Rusin 2011), and the possibility of non-solar abundances 
must be considered (Reynolds, Fabian & Inoue 1995). 
While acknowledging our uncertainties about the nature 
of the primary X-ray source in AGN, the physicality of 
the irradiation profile e(r) across the disk must be re- 
spected. For example, extremely centrally concentrated 
irradiation profiles (e oc r~ q ,q > 5) are difficult to un- 
derstand unless the black hole is rapidly-rotating. Spec- 
tral variability presents yet another complicating factor, 
often putting us in the position of having to model a 
time-averaged spectrum that has been collected over a 
time-frame where we know that spectral variability has 
occurred. Finally, of course, attention must be paid 
to modeling the photoelectric absorption features that 
sometimes dominate the soft X-ray spectra of AGN and 
may itself be time- variable. 

In this paper, we revisit the case of the bright Seyfert 
1.5 galaxy NGC 3783 (z = 0.00973; Theureau et al. 
1998). Using a deep (250 ks) Suzaku observation, Bren- 
neman et al. (2011; hereafter Bll) identified strongly 
broadened spectral features and concluded that the black 
hole is rapidly-spinning in the prograde sense (with re- 
spect to the accretion disk), with a dimensionless spin 
parameter a > 0.88 at the 99% confidence level. Reis et 
al. (2012) examined spectral variability during this same 
observation. They found that, while the strong warm ab- 
sorber showed no signs of variability, there were signifi- 
cant intra-observation changes in the 15-40keV/2-10keV 
color which they modeled as ionization changes within 
the disk. However, they showed that the overall con- 
clusions of Bll were unaffected by this variability and, 
again, a rapidly rotating black hole was required. How- 
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ever, using the same dataset, Patrick et al. (2011; here- 
after Pll) obtained a different result, preferring a slowly 
or retrograde spinning black hole a < —0.04. Clearly, 
this discrepancy must be addressed and the robustness 
of the spin in this (and other AGN) assessed. 

We identify the treatment of the iron abundance as 
the key that underlies the discrepancies in the published 
analyses of NGC 3783; Pll assume that the accretion 
disk possesses cosmic iron abundance whereas Bll allow 
for non-solar abundances. Revisiting these data with a 
Monte Carlo Markov Chain (MCMC) analysis, we find 
that black hole spin and iron abundance are statistically 
correlated variables. We show that the data strongly sup- 
port a super-solar iron abundance (Z-p = 2.3 — A.IZq; 
90% confidence) and a rapidly spinning black hole (a > 
0.89; 90% confidence) in this object. While the iron 
abundance found in NGC 3783 is completely in line with 
the general metallicity enhancements of AGN inferred 
from optical/UV spectroscopy, we note that the X-ray 
reflection technique has found extreme iron abundances 
in several narrow line Seyfert 1 galaxies. We suggest that 
radiative levitation can enrich the photosphere of the in- 
ner accretion disk with iron, leading to apparent extreme 
iron abundances. 

This paper is organized as follows. Section 2 gives 
a brief description of the Suzaku data and its reduc- 
tion. Section 3 describes the spectral modeling of these 
data, and the results from this modeling are presented 
in Section 4. The implications of our study are dis- 
cussed Section 5, where we also make a preliminary as- 
sessment of the ability of radiative levitation to enhance 
the photospheric abundances of the innermost accretion 
disk. Throughout this paper, luminosities and distances 
are calculated using a ACDM cosmological model with 
H = 71kms~ 1 Mpc" 1 , fi A = 0.73 and Q M = 0.27 (Ko- 
matsu et al. 2011). For a redshift of z = 0.00973, 
this results in a luminosity distance to NGC 3783 of 
41.4 Mpc. Unless otherwise stated, all errors in this pa- 
per are quoted at the 90% confidence level for one inter- 
esting parameter. 

2. OBSERVATION AND DATA REDUCTION 

Suzaku observed NGC3783 quasi-continuously for the 
period 2009 July 10-15. The source was placed at the 
nominal aimpoint for the Hard X-ray Detector (HXD). 
After filtering the data for Earth occultations, South 
Atlantic Anomaly passages, and other high background 
events, we obtained a total on-source exposure of 210 ks 
spread over 360 ks of wall-clock time. The data reduc- 
tion exactly follows that described in Bll except for the 
use of updated calibration and background files (as of 
late-Fall 2011). We use data from all three operating 
X-ray Imaging Spectrometers (XIS0/XIS1/XIS3), com- 
bining the spectra from the two front illuminated (FI) 
detectors (XIS0 and XIS3) into a single spectrum. When 
fitting spectra, the flux cross-normalization between the 
XIS-FI spectrum and the back illuminated (XIS1) spec- 
trum is left as a free parameter and is found to be ap- 
proximately 1.03 (the same as with the older calibra- 
tions used by Bll). Hard X-ray coverage is provided 
by the HXD/PIN instrument, reduced as per the Suzaku 
ABC guide using the relevant "tuned" non X-ray back- 
ground event file. During spectral fitting, the PIN/XIS- 
FI flux cross-normalization was fixed at 1.18 (Suzaku 



Memo 2008-063). In Bll, we investigate the effects of 
allowing this cross-normalization to be a free parameter, 
finding that it adopts a value of 1.15 ± 0.07 with rather 
little impact on other aspects of the fit. 

The spectral analyses presented in this paper are per- 
formed jointly on the simultaneous XIS-FI, XIS1, and 
HXD/PIN spectra. We use XIS data in the 0.7-10 keV 
band, excluding the 1.5-2.5 keV region due to the pres- 
ence of well known calibration artifacts associated with 
the X-ray mirror and Si-edge in the detector. XIS data 
were rebinned to a mininum of 20 photons per energy bin 
to facilitate the use of % 2 statistics, and PIN data were 
regrouped to a minimum signal-to- noise ratio of 5. 

3. SPECTRAL MODELING 

The spectrum of NGC 3783 is complex and its descrip- 
tion requires a multi-component model. We begin by 
discussing the basic components and parameters of the 
base model. At the end of this section, we delineate the 
specific variants of the model that we shall employ. 

The primary X-ray continuum is described by a power- 
law with photon index T. X-ray reflection from distant 
neutral material (the outer accretion disk or the dust 
molecular torus) is described using the pexmon model 
(Nandra et al. 2007; Bll) that includes the iron Ka/K/3 
fluorescent emission lines and the associated Compton 
reflection continuum. The strength of this neutral reflec- 
tion is characterized by the reflection fraction, 1Z, nor- 
malized such that 1Z — 1 corresponds to a reflector that 
subtends a solid angle of 27T as seen from the source. 
To model an additional narrow feature in the iron K- 
band, we need to add a narrow Fe XXVI emission line 
at 6.97 keV which presumably arises from highly pho- 
toionized optically-thin matter close to the AGN. X-ray 
reflection from the ionized inner accretion disk is de- 
scribed with the ref lionx model (Ross & Fabian 2005) 
convolved with the relconv kernel in order to include 
relativistic Doppler and gravitational redshifts (Dauser 
et al. 2010). The rest-frame reflection spectrum of the 
disk is characterized by its ionization parameter £, and 
the iron abundance relative to cosmic Zp^disk- The rela- 
tivistic smearing of the disk spectrum depends upon the 
spin of the black hole a, the inclination of the inner disk 
i, the radius of the inner edge of the accretion disk (set to 
the innermost stable circular orbit [ISCO] for that spin), 
and the irradiation/emissivity profile of the disk (set to 
be a broken power law, breaking from r~ qi to r~ q2 at 
radius r = rb r ). Provided that (72 > 2, the line profile 
is insensitive to a (sufficiently large) outer radius; we fix 
fout = 400r ff , where r g = GM/c 2 is the gravitational 
radius. 

The multi-component spectrum of NGC3783 is 
strongly affected by photoionized absorption which, 
based upon Chandra/ HETG data (Kaspi et al. 2002, 
Netzer et al. 2003, Brenneman et al. 2011), we model 
using a three-zone warm absorber (with column densi- 
ties Ni and £j, i — 1,2,3). We allow a fraction / of 
the continuum emission to scatter/leak around the pho- 
toionized absorbers. We use the XSTAR code (Kallman 
& Bautista 2001) to compute a grid of absorption spectra 
on the 2-dimensional (£, iV)-plane, sampling the param- 

7 http: / / www. astro. isas . j axa.jp/suzaku /doc /suzakumemo /suzakumemo- 
2008-06.pdf 



MCMC investigation of spin in NGC3783 



3 




Figure 1. Illustration of the importance of £- resolution in con- 
structing warm absorber tables. The solid lines show the abun- 
dance of four important ions (O VII, O VIII, Si XIII, and Sixiv) rel- 
ative to the total abundance of that element as a function of the 
ionization parameter £; these curves are computing using XSTAR 
assuming a plasma density of 10 10 cm -3 and a power-law ioniz- 
ing spectrum with photon index T = 2 extending from 1 Ry to 
1000 Ry. The squares connected by dotted lines illustrate the ac- 
curacy with which linear interpolation can reproduce these curves 
with a 5-samples-per-decade criterion (A log J = 0.2). The dashed 
lines illustrate a linear interpolation using just 2. 5-samples-per- 
decade. The lower resolution sampling leads to substantial errors 
in the reconstructed ion abundance (and hence opacities). 

eter space at 20 logarithmically-spaced points in column 
density (from N = 10 20 cnr 2 to N = 10 24 cm~ 2 ) and 
20 logarithmically-space points in ionization parameter 
(from £ = lergcms -1 to £ = 10 4 ergems -1 ). For the 
purpose of defining the ionization balance in the absorb- 
ing zone, the photoionizing spectrum is fixed at T = 2, 
a good approximation to the actual primary X-ray spec- 
trum of NGC3783. The velocity dispersion of the ab- 
sorbing gas is set to 200 kms -1 . Motivated by Pll, 
we experimented with using a higher velocity dispersion 
(1000 kins" 1 ) for the most highly ionized absorber but 
found that the results were unchanged. 

We note that, when employing tabulated models to 
describe a strong warm absorber, it is crucial to sample 
the parameter space with sufficient resolution. Of spe- 
cial importance is the ^-resolution of the tabulated grid. 
The abundance of a given ion (say O vn or Si XIV) dis- 
plays a broad peak as a function of £. This peak must 
be resolved by the tabulated grid if one wishes to obtain 
accurate spectra when interpolating to intermediate val- 
ues of £. This is illustrated in Fig. [TJ which shows the 
ion abundance (relative to the total abundance of that 
element) for Ovu, Ovin, Sixm, and Sixiv as a func- 
tion of ionization parameter £, calculated using XSTAR 
assuming a plasma density of 10 10 cm -3 and a power- 
law ionizing spectrum with photon index T — 2 extend- 
ing from 1 Ry to 1000 Ry. Sampling £ with a resolution 
Alog£ = 0.2 and linearly extrapolating between sam- 
ple points (as does XSPEC) produces estimates of the 
ion abundances accurate to better than 5%, translating 



directly into interpolated warm absorber grids that will 
possess opacities accurate to better than 5%. These con- 
siderations drive the construction of our absorber grids; 
we have 5 models per decade of £ across 4 decades of 
£. It does raise serious concerns about results obtained 
using the XSPEC model zxipcf; the grid underlying this 
model only samples 12 points across 9 decades of £ and 
hence is clearly undersampled. This resolution issue may 
also be another difference between Bll and Pll; the XS- 
TAR tables used for the Pll analysis either sample 5 
decades of £ with 11 values (for the low- and moderate- 
ionization absorbers) or 6 decades of £ with 20 values (for 
the high- ionization state absorber). Thus, the sampling 
falls somewhat below our 5-samples-per-decade criterion 
and interpolation to intermediate ionization states may 
become inaccurate. 

Finally, we find the need to include a soft excess com- 
ponent — following previous works, we employ a purely 
phenomenological approach. Previous treatments of the 
soft excess in NGC 3783 have used a blackbody compo- 
nent (characterized by temperature Xbt,) or Comptoniza- 
tion by "warm" plasma (described using the comptt 
model in XSPECV12.7.0 and characterized by a tempera- 
ture T comp t and an optical depth r). To be conservative, 
we initially included both a blackbody and a warm ther- 
mal Comptonization component in our spectral fit. In all 
cases presented in this paper, the Comptonization com- 
ponent would drop out of the fit (with its normalization 
becoming zero). Experimenting with including either 
the blackbody or the Comptonization component suggest 
that the blackbody fit was slightly preferred (A% 2 ~ 3); 
in both cases, the components have essentially the same 
shape and contribute only to the lowest energies in the 
Suzaku energy band. Henceforth, all of the spectral mod- 
eling in this paper uses a blackbody model for the soft 
excess. Using spot-checks, we have confirmed that the 
key results (e.g. the conclusions regarding black hole 
spin) are unchanged by switching to a Comptonization 
realization for the soft excess. 

Here, we shall employ three realizations of this spec- 
tral model that differ in their treatment of iron abun- 
dance. In Model A, we assume that the circumnuclear 
region (inner accretion disk, outer accretion disk, and 
the obscuring torus) are chemically homogeneous. We 
make no assumption, however, about the value of this 
common iron abundance; we allow it to be a free pa- 
rameter (Zpodisk = ^Fctor = free). Model B also as- 
sumes a chemically homogeneous system but imposes the 
additional constraint that the iron abundance is solar 
(Zpe.disk = -^Fe.tor = 1)- Model C describes a chemically 
inhomogeneous system in which the distant reflector (the 
torus or the outer accretion disk) are forced to have so- 
lar iron abundance but the inner accretion disk has a free 
iron abundance (^Fe^isk = free, ^Fc,tor = 1). 

4. RESULTS 

4.1. The need for super-solar iron abundances and rapid 
black hole spin 

Model A is essentially identical to the spectral model of 
Bll and a traditional y 2 analysis of the Suzaku spectrum 
with this model was published in that paper. As already 
noted in the introduction, Bll finds a rapidly spinning 
black hole and super-solar iron abundance. However, a 
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Figure 2. MCMC results for black hole spin and iron abundance from Model A. Top panels : Parameter value along the 250,000-element 
chain trajectory for the black hole spin a (left) and iron abundance Zp e (right). For clarity, we only show every 100th element of the chain. 
Bottom panels : Probability distributions for the black hole spin and iron abundance. 

significant concern when x 2 fitting a spectral model with 
a large number of free parameters (N = 22 free parame- 
ters in our case, including all model normalizations and 
the FI/XIS1 cross-normalization factor) is that we will 
mistake a local minimum in x 2 space for the global best 
fit. Furthermore, the spectral model may genuinely pos- 
sess multiple solutions and traditional x 2 error analysis is 
not well suited for mapping out multi-peaked probability 
distributions for model parameters. 

To address these limitations and give an independent 
validation of the statistical analysis of Bll, we employ 
MCMC techniques which are now implemented in the 
Xspec spectral fitting package. Operationally, we pro- 
ceed as follows. Using traditional x 2 minimization, we 
obtain the best fit of Model A with the full multi- 
instrument Suzaku dataset. We then run five 55,000 ele- 
ment chains, starting from a random perturbation away 
from the best- fit and ignoring ( "burning" ) the first 5000 
elements of each chain. The proposal for the chain (i.e. 
probability distribution that dictates the Monte Carlo 



step for each model parameter as the chain proceeds) is 
drawn from the diagonal of the covariance matrix result- 
ing from the initial x 2 minimization, assuming Gaussian 
probability distributions with a rescaling factor of 10~ 3 . 

Figure |5] (top panels) shows the values of the spin 
and iron abundance through the 250,000 elements of the 
chains Model A, showing every 100th element for clar- 
ity. The lack of large excursions or trends suggests that 
we have indeed found the global best fit of this spec- 
tral model. Experiments in which we start at a distant 
point in parameter space (e.g., setting initial seeds of 
a = —0.95 and Zp e = 1) converge on the same best fit. 
Similarly, employing a different proposal (using Cauchy 
probability distribution or a larger rescaling factor) also 
converges on the same best fit, albeit slowly and with a 
very high repeatability fraction in the chain. 

The best-fit and MCMC-derived (90% confidence level) 
errors are reported in Table [TJ and Fig. [2] (bottom pan- 
els) reports the detailed probability distributions for 
the black hole spin and iron abundance resulting from 
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Table 1 

Spectral parameters for the model fits discussed in this paper. Column densities are in units of 10 20 cm -2 . Ionization parameters are in 

cgs units. The power-law normalization is in units of 10 -2 phs -1 keV measured at 1 keV. The blackbody normalization is in units of 
10 phs -1 , and its temperature is in units of keV. The energy and 1-a width of the Fe26 line is in units of keV and the equivalent width 
Wfe26 is in units of eV. Fixed parameters are denoted by "(f)". All errors are quoted at the 90% confidence level. 



Model A. Results are, within the error bars, completely 
consistent with the results of Blfl From the relativisti- 
cally smeared disk features, we derive a black hole spin 
of a > 0.89, disk inclination between i = 21 — 26°, and 
an iron abundance between Z — 2.3 — A.IZq. The in- 
ner emissivity index is steep (gi — 6.6^5 '4), in agreement 
with Bll and consistent with a rapidly spinning black 
hole. The emissivity profile flattens to <j2 = 2.9 ± 0.2 be- 
yond the break radius (r^ r = 3.7 — 5.5r g ). Thus, without 
any prior constraint or bias, the outer emissivity index 
is consistent with the fiducial value of qi = 3 as pre- 
dicted by either a lamp-post illumination model or any 
model in which the X-ray irradiation/emissivity tracks 
the underlying dissipation in the accretion disk. 

An advantage of MCMC techniques is the efficiency 
with which statistical correlations between parameters 
can be uncovered. Here, we are particularly interested 
in parameters that statistically correlate with black hole 

8 We note that goodness of fit for this model is formally un- 
acceptable (x 2 /dof = 1.33 for 687 de grees of freedom). However, 
detailed visual inspection of the model fit to the folded spectrum 
shows that this is due to residuals in the 0.7-1.5 keV range; re- 
moving these energies from consideration results in a formally ac- 
ceptable fit x 2 /dof = 1.03 (556 degrees of freedom). The fact that 
the soft X-ray residual features are narrower in energy than the 
instrumental response, and are uncorrelated between XIS-FI and 
XIS1, strongly points to calibration features. However, the small 
amplitude of these residuals (< 5%) and the fact that the model 
averages between the two XIS spectra gives us confidence that they 
are not dramatically skewing our best fit model. 



spin. The strongest such correlation is between spin and 
iron abundance. We illustrate this coupling in Fig. [3] 
(left) which shows the 2-dimensional probability distri- 
bution on the (Zpc, a)-plane for our Model A analysis. 
There is clearly a positive correlation between these two 
parameters, with more rapid spins simultaneously favor- 
ing higher iron abundance. It is straightforward to un- 
derstand this correlation. As we discuss further in Sec- 
tion 231 the sensitivity to spin comes as much from the 
need not to over-predict the blue horn of the iron line as 
it does from the need to capture the red-wing of the line. 
As the iron abundance of the model increases, the ratio 
of iron line photons to reflection continuum increases. 
In turn, the fit must blur the iron line more severely in 
order to not over- model the data in the 5-6 keV band, 
and hence the inferred spin increases. Similar reasoning 
explains the weak anti-correlation between the reflection 
fraction of the distant reflector 1Z and the black hole spin 
(Fig. El middle). As 1Z decreases, the contribution of 
the inner accretion disk to the observed reflection spec- 
trum increases leading to more iron line photons from 
the disk. Consequently, the model needs to blur the disk 
line more severely (pushing to higher spin) in order to 
not over-model the data in the 5-6 keV band. 

A less obvious statistical coupling is between the spin 
and and the ionization state of the intermediate warm 
absorber (£2; Fig. [3] right). The 90% confidence range 
of this parameter is log £2 = 1-85 — 1.92 (cgs units). We 
find that the spread in allowable black hole spins is sig- 
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Figure 3. 2-dimcnsional probability density functions of those parameters that show a correlation with spin. Left panel : Probability 
density on the (Zp e , a)-plane showing a positive correlation. Contour levels are shown at p(Zp e ,a) = 1,3.3,10.0 [defined such that the 
probability of being in the range Zp G — ► Zp G + <5Zp c and a — > a + Sa is p(Zp c , a)SZp c Sa]. Middle panel : Probability density on the 
a)-plane revealing a weak negative correlation. Contour levels are shown at p(1Z, a) = 10, 33, 100. Right panel : Probability density on 
the (log £, a)-plane revealing a weak negative correlation. Contour levels are shown at p(log £, a) = 10, 33, 100. 



nificantly greater when log £2 > 1-9, resulting in a "T" 
shape to the 2-d probability on the (log £2, a)-plane. No 
similar coupling is seen with £1 or £3. We conclude that, 
in the range log £2 = 1.9 — 2.0, the spectral curvature 
can be (in part) degenerately described by either strong 
relativistic blurring characteristic of a very rapidly rotat- 
ing black hole or by absorption edges from this the warm 
absorber (see Young et al. 2005 for a detailed discus- 
sion of these "broad line mimicking warm absorbers"). 
The important thing to note is that, even accounting for 
all of these statistical correlations and degeneracies, the 
spectral fit still strongly requires a rapidly rotating black 
hole. 

In addition to these statistical correlations, we note 
that both the 1-d and 2-d probability distributions 
(Figs. [2] and [3]) show double maxima, with peaks at 
a fn 0.93 and a — 0.99. We have confirmed that the 
relativistic convolution model varies smoothly in the ex- 
pected manner across this parameter range and, hence, 
is not an obvious source for this bi-modality. We sus- 
pect that small (< 1%) systematic errors in the effective 
area models (especially at energies characteristic of the 
iron line red-wing, 3-4 keV) have induced this structure 
in probability space. 

4.2. Models with imposed abundance constraints 

Having noted the existence of a correlation between 
iron abundance and spin, we now make a connection 
with the work of Pll and explore a model (Model B) in 
which solar abundances are imposed. We have repeated 
our MCMC analysis, now employing Model B, again run- 
ning five 55,000 element chains with a 5000 clement burn 
period. 

Imposing this abundance constraint (and refitting all 
other parameters) worsens the goodness of fit parameter 
by Ax 2 — 36 (for one fewer degree of freedom) and sig- 
nificantly changes the probability distribution for several 
of the key spectral parameters (Table [1} . The principle 
effect is that, due to the decreased iron abundance, the 
fit demands a significantly higher neutral reflection frac- 
tion 1Z in order to fit the narrow iron line. This results 
in a slightly steeper best-fitting power-law continuum (r 
steepens from 1.79 to 1.85). At the same time, the ratio 
of iron line photons to reflection continuum from the ac- 



cretion disk is also decreased by the reduced abundance. 
Thus, the need to severely broaden the iron line and, 
hence, invoke a rapidly spinning black hole, is dimin- 
ished. In detail, the probability distribution for spin in 
this constrained fit is weakly bi-modal, with modest-to- 
rapid prograde or retrograde black holes being preferred 
over slowly spinning holes (Fig. |4]). 

Both Model A and Model B assume that the distant 
reflector and the inner accretion disk have the same iron 
abundance. In our final spectral model, Model C, we 
explore the effects of imposing a mixed abundance model 
such that the distant reflector has a solar iron abundance 
but we allow the inner accretion disk to have a freely- 
fitting iron abundance. In terms of goodness of fit, this is 
intermediate between Model A and Model B, with A% 2 = 
14 compared with our fiducial fit. The MCMC analysis 
(executed in an identical manner to that for Model A) 
shows that a rapidly spinning black hole is preferred, but 
there is a low-spin tail to the probability distribution 
giving a final result of a > 0.60 (at 90% CL). 

4.3. What drives the spin constraint? 

Fundamentally, analyses of disk reflection features di- 
agnose black hole spin using the fact that the observed 
broadening/redshifting can be so extreme as to demand 
reflection from within 6r g , the ISCO of a non-rotating 
black hole (Reynolds & Fabian 2008). Thus, one might 
expect that the spin constraints come principally from 
the modeling of the extreme red-wing of the iron line 
and that, if the spectral model being compared with the 
data has too low of a spin, the main result would be pos- 
itive residuals (model under-predicting data) at energies 
characteristic of the red-wing (3-5 keV). In fact the sit- 
uation is more subtle due to the ability of other model 
components and parameters (e.g. the warm absorber, or 
the iron abundance of the accretion disk) to compensate 
for an incorrect spin. 

We illustrate this in Fig.[5J Figure^ shows the model 
and residuals for the best fit model presented in Sec- 
tion 14.11 For comparison, Fig. [5}d shows the model and 
residuals for a fit in which the spin has been fixed to a = 
and the emissivity indices have been fixed to q\ — (72 = 3 
(to be physically consistent with a non-spinning black 
hole) . While there are weak positive residuals in extreme 
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Figure 4. Probability distribution for black hole spin in our fidu- 
cial model (Model A, blue), the solar abundance model (Model B, 
red), and the inhomogeneous-abundance model (Model C, green). 

wing of the iron line at 3-4 keV, the most significant de- 
viations are negative residuals at 5-6 keV. This can be 
understood as follows. The overall strength of the ac- 
cretion disk reflection as well as the iron abundance of 
the disk is, to a major extent, determined by the iron- 
edge (rest frame energy 7.1 keV) and the form of the hard 
band (7-40 keV) spectrum. The question then becomes 
whether an iron line with the self-consistently determined 
equivalent width is compatible with the data. In the 
case of NGC3783, we find the resulting iron line is too 
strong (over-predicts the data) if it is subjected to blur- 
ring from around a non-spinning black hole. Only the 
enhanced relativistic effects of a rapidly-rotating black 
hole can distribute the line flux sufficiently to make the 
iron line compatible with the data. 

Figure[5]3 teaches us an important lesson about the role 
of possible systematic errors. Specifically, given data of 
realistic complexity, the relative effective area calibration 
of a medium-resolution instrument such as the XIS must 
be accurate to better than 2% across the 3-7 keV band 
if we hope to detect the signatures of black hole spin. A 
more complete study of the role of systematic errors in 
black hole spin measurements (including an assessment 
of how the effective area calibration requirements relax 
when one has simultaneous high-resolution constraints 
on the warm absorber, as will be the case for the Soft 
X-ray Spectrometer [SXS] on Astro-H) will be presented 
in a future publication (A.Lohfink in prep.). 

5. DISCUSSION 

5.1. Summary of observational results 

We have revisited the long 2009 Suzaku data for the 
Seyfert 1.5 galaxy NGC 3783 in the light of conflict- 
ing published results for the black hole spin. We find 
that the treatment of the iron abundances in both the 
inner accretion disk and the distant torus are impor- 
tant for the derived SMBH spin. We expect that the 
timescale on which stellar processes enrich their environ- 



ment (~ 10 s yr) is longer than the timescale on which 
material is exchanged between the torus and the in- 
ner disk (of order the dynamical timescale of the torus, 
10 5 yr) leading to the conclusion that the central engine 
will be chemically homogeneous. Employing this ho- 
mogeneity assumption (Zpe.disk = ^Fc,tor), our spectral 
analysis shows that the data strongly require high iron 
abundance Zp c = 2 — 4Z Q and a rapidly spinning black 
hole a > 0.88. We note that such a metallicity enhance- 
ment is completely in line with that inferred from opti- 
cal/UV studies of the broad line region (BLR) of type-1 
AGN, e.g. see Warner et al. (2004) and Nagao, Marconi 
& Maiolino (2006). In agreement with previous analy- 
ses, we show that constraining the abundance to be solar 
worsens the goodness of fit and leaves the SMBH spin 
unconstrained. Employing a chemically-inhomogeneous 
model (Zpctor = Z Q ; ZF Gj( jisk =free) we show that it 
is the abundance of the accretion disk component that 
is primarily (but not exclusively) driving the spin con- 
straint. 

5.2. Iron abundance measurements from X-ray 
reflection 

Even if the accretion disk and the surrounding circum- 
nuclear region are chemically-homogeneous, there are 
subtleties that can affect the abundances that are in- 
ferred from X-ray reflection fits. Firstly let's consider 
the "cold" reflection from the torus. Fundamentally, X- 
ray reflection models constrain iron abundance by com- 
paring the equivalent width of the fluorescent iron line 
to the strength of the Compton reflection continuum at 
E > 7keV. With photons at these energies, the nature 
of the Compton scattering and hence the reflection con- 
tinuum is the same whether the bulk of the electrons are 
free or bound into hydrogen/helium atoms. However, of 
potential importance is the depletion of iron into dust 
grains within the torus. While an iron atom can read- 
ily undergo K-shell fluorescence even when bound into a 
solid, sufficiently large dust grains will lead to the self- 
shielding of iron and hence a reduction in the iron-Ka 
equivalent width. The minimum grain size a; m i n for which 
K-shell self-shielding is important is given by assuming 
pure iron, x m i n = TOFe/pFeCth where mpe is the mass of a 
iron atom, pp c is the density of the grain, and o~th is the 
cross-section at the K-shell photoionization threshold for 
cold iron. Evaluating, using ppo = 7.7 x 10 3 kgm -3 and 
o~th — 3.5 x 10~ 24 m -2 , we obtain x m - m k 3.5 /iin. The 
minimum grain size may be up to twice this size is the 
iron if bound into a compound (e.g. an oxide, silicate or 
sulphide), and even more if the grain has an amorphous 
structure. Thus, self-shielding can only reduce the iron- 
Ka equivalent width if the dust mass is dominated by 
very large grains. 

For the innermost accretion disk, different considera- 
tions are at play. Since there are clearly no issues of 
dust depletion in this hot plasma, the normal assump- 
tion is that the abundances measured for the inner ac- 
cretion disk should faithfully represent the true abun- 
dances. However, in this extreme environment, there are 
some interesting mechanisms that have the potential to 
affect the photospheric iron abundance of the accretion 
disk. Skibo (1997; also see Turner & Miller 2010) sug- 
gests that the bombardment of the inner accretion disk 
by energetic (relativistic) protons can lead to the spalla- 
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Figure 5. Left panel : Unfolded XIS data, best fitting model and resulting featureless residuals for Model A. The model components 
visible here are the absorbed power-law continuum (red), distant reflection (green) and relativistically smeared disk reflection (blue). Right 
panel : Unfolded XIS data, best fitting model and residuals when the spin is constrained to be a = and the emissivity indices are 
constrained to be q± = q2 = 3 (but all other parameters are refitted). Model components are colored the same as above. 



tion of iron, enhancing the abundance of Ti, V, Cr and 
Mn at the expense of the iron abundance. However, sig- 
nificant iron spallation within the disk only occurs if a 
large fraction of the accretion energy (~ O.lMc 2 ) is car- 
ried by relativistic protons which subsequently bombard 
the disk. The resulting 7-ray emission (resulting from 
the decay of pions produced by these interactions) was 
estimated by Skibo (1997) and violates the 4cr upper lim- 
its for radio-quiet Seyfert galaxies recently found Teng et 
al. (2011) using Fermi data. 

Most AGN in which X-ray reflection from the inner ac- 
cretion disk has been modeled in detail have been found 
to possess supersolar iron abundance (e.g., see Table 2 
in Bll). In most cases, such as NGC 3783, the inferred 
abundances are completely consistent with that found for 
the BLR and hence probably reflects a true-chemically 
enriched circumnuclear environment. In a small number 
of cases, however, rather extreme disk iron abundances 
are inferred (e.g., 1H0707— 495 in which a strong iron- 
L line implies iron abundances of > 7Z@; Zoghbi et al. 
2010). In the rest of this Section, we discuss a mecha- 
nism, radiative levitation, that has the potential to sig- 
nificantly enhance the photospheric iron abundance of 
the inner disk. 

5.3. Radiative levitation in AGN accretion disks 

In moderate luminosity AGN such as NGC 3783, 
the inner regions of the accretion disk are radiation- 
pressure dominated with a photospheric temperature of 
T ~ 10 5 - 10 6 K (Shakura & Sunyaev 1973). By defini- 
tion, this implies that each radius of the disk is at the 
"local" Eddington limit in the sense that the radiation 
pressure force acting on the plasma is balancing the ver- 
tical component of gravity. However, at the microphys- 
ical level, the radiative force on metal ions can be very 
different to the radiative force on the background, fully 
ionized, hydrogen/helium plasma. In particular, since 
iron can possess moderate-to-low charge states (Fe XVII 



and below) with populated L- and M-shells, the radia- 
tion force of the 10 5 — 10 6 K thermal radiation pushing 
on the UV/EUV iron absorption lines and the EUV/X- 
ray bound-free edges can produce a net upwards force on 
the iron ions many orders of magnitude greater than the 
vertical gravity. In the absence of other dynamics (but 
see below), this radiative pushing will cause the iron to 
diffuse upwards towards the photosphere and could sig- 
nificantly enhance the iron abundance of the disk sur- 
face. This phenomenon has been discussed in the context 
of the surface abundances of hot white dwarfs (Chayer, 
Fontaine & Wesemael 1995; Seaton 1996, Wassermann 
et al. 2010), and opacity enhancement from radiative 
levitation have been identified as key to understanding 
the 1 2 min oscillations in the unusual sub-dwarf O-star 
SDSS J160043.6+074802.9 (Fontaine et al. 2008). How- 
ever, to the best of our knowledge, radiative levitation 
has not been discussed previously in the context of AGN 
accretion disks. Here, we conduct a preliminary assess- 
ment of whether radiative levitation can be relevant for 
AGN accretion disks — we conclude that it may be rel- 
evant, and lay out a path for a more rigorous treatment 
of its effect. 

Consider an iron ion (mass mpe and charge state Qp c ) 
that is within a few Thomson depths of the photosphere 
of the inner accretion disk at radius R = xr g from the 
black hole. The dimensionless accretion rate of the ac- 
cretion disk is taken to be to = L/rjL^d where L^dd 
is the Eddington luminosity and rj is the accretion ef- 
ficiency. The standard theory of radiation-dominated 
disks (Shakura & Sunyaev 1973) then gives the disk 
thickness as h = ffTOV Our iron ion is subject to a 
gravitational force pulling it towards the midplane of the 
disk, -F grav = TOFeS where g — GMh/R 3 . The ion is also 
subject to a radiative force pushing it upwards towards 
the photosphere F ra d = TOF C g ra dj as well as the effect 
of Coulomb collisions with the ions in the background 
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plasma (assumed, for simplicity, to be pure hydrogen). 
Following Seaton (1996), the (upward) drift velocity of 
the iron ion relative to the background plasma is given 

by 



V 



(ffrad - g) 



ldx 



(1) 



where T is the temperature of the plasma (assumed to 
be in LTE), \ is the abundance enhancement factor of 
the iron ions, and V the diffusion coefficient. For this 
non-degenerate plasma, the appropriate form of the dif- 
fusion coefficient is (Chapman & Cowling 1962; Alcock 
& Illarionov 1980), 



V 



3(2fcT) 5 / 2 
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where n is the proton number density of the background 
hydrogen plasma and A is the Coulomb logarithm given 

by, 

(kTf 
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To assess the potential relevance of radiative levita- 
tion, we will estimate the drift timescale across a Thom- 
son depth, Adrift = 1/tuttV. Competing mechanisms (e.g. 
"stirring" of the surface layers) will occur on the order 

of the dynamical timescale td yn — \J R 3 /GM, and thus 
we shall examine the ratio of timescales idriftAdyn- We 
will also limit our consideration to the situation where 
dx/dz — (i.e. no enhancement has yet occurred). Us- 
ing the above expressions and definitions, a little algebra 
shows that 



tdrift _ 16 (vrm^/VA Q 



,3/2 



i dyn 9 m Fc a T c (2fcT) 3 / 2 mS rad : 



where we have defined, 
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Note that, because we are considering the drift timescale 
across a Thomson depth, the proton density of the 
plasma n only appears in the expression for tdriftAdyn 
via the Coulomb logarithm A; hence the result is ex- 
tremely weakly dependent on the density. Now, suppose 
that the net upward energy flux of radiation (as a func- 
tion of frequency v) is F v . The net radiation force on the 
iron ion is 

1 f°° 

™Fc.9rad = - / O v F v dv, (6) 

c Jo 

where a v is the absorption cross-section of the ion. As- 
suming that the opacity of the disk is dominated by elec- 
tron scattering, we can use the fact that the disk is radi- 
ation pressure supported to write, 



<7 f°° 

m p g = — / F v dv, 
c Jo 

which, when combined with eqn. |6] gives, 
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Figure 6. Estimates of the radiative levitation drift timescale for 
various charge states of iron (expressed as the ratio tdriftAdyn X 
mi" 3 ' 2 ), calculated under the assumptions laid out in £|5.3I Shown 
here is Feviii (black; solid), Fbxvii (red; solid), Fexviii (red; 
dotted), Fexix (red; dashed), Fexx (red; dot-dashed), and Fexxi 
(red; dot-dot-dot-dashed). 

A full evaluation of F v and hence S ra d must be per- 
formed in the context of a detailed radiative transfer 
model; such a treatment is beyond the scope of this pa- 
per. Here, we make an enormous simplifying assump- 
tion, that the spectral shape of F v is proportional to 
the Planck function £?„(T). While this is certainly in- 
appropriate, it does allow us to study whether radiative 
levitation has the potential to be important. We also as- 
sume that only opacities from bound-bound transitions 
are relevant. Under these assumptions, a little algebra 
shows that 



' — rad 



45 m e m p c h 
87r 4 mF C e 2 kT 



fi 



Hi 



1 



(9) 



where we have replaced the integral in eqn. [6] with 
the sum over bound-bound transitions with oscillator 
strengths fa and normalized energies yi = hvi/kT . Given 
the appropriate atomic data, we now have all of the in- 
gredients required to compute idriftAdyn as a function 
of charge-state Qfo, temperature T, dimcnsionless accre- 
tion rate rh and dimcnsionless radius x. 

Figure [S] shows this ratio tdriftAdyn (normalized to 
x = m — 1) as a function of temperature for vari- 
ous charge states of iron. Specifically, we show the re- 
sults for Feviii (the dominant charge state of iron for 
log£ = 0.4 - 1.2) and Fexvii-Fexxi (the dominant 
charge states for log£ = 1.8 — 2.8). The atomic data 
needed to evaluate eqn. [S] is obtained from the the XS- 
TAR database. It can be seen that, for T > 3 x 10 5 K, 
the drift timescale for Feviii and Fexvii can be only 
a factor of a few greater than the dynamical timescale. 
For Fe xviii-Fe XXI ions, the drift timescale is approxi- 
mately an order of magnitude greater than the dynamical 
timescale for T = 10 6 K, but becomes rather longer for 
smaller temperatures. We conclude that radiative levi- 
tation may be an important process in the inner disks 
(x ~ 1 — few) of high accretion rate (m ~ 1) low black 
hole mass (high disk temperature) AGN provided that 
the ionization state promotes certain charge states of 
iron. 
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The above calculation is only a plausibility argument. 
Firstly, as mentioned already, a real radiative transfer 
calculation is required to compute the radiative force on 
the iron ions. Line saturation will be extremely impor- 
tant and tend to decrease the upward radiative force, but 
this may be partially counteracted by ionization gradi- 
ents and Doppler shifts associated with the radiation in- 
duced ion drift. We also note that our calculation does 
not include the radiation force resulting from pushing 
on bound-free transitions for which saturation will be 
much less important. Given sufficient time (many drift 
timescales) , radiative levitation would achieve an equilib- 
rium in which the iron enhancement produces increased 
opacity (self-shielding) such as to bring the radiative and 
gravitational forces into balance (Chayer, Fontane & We- 
semael 1995). But, as a second major issue, we must 
characterize processes that will mix the plasma and erase 
the radiation- induced enhancements. If the surface layers 
are turbulent with eddy timescales comparable to or less 
than the local dynamical timescale, we see that radiative 
levitation may be of limited importance. On the other 
hand, if the surface layers of the disk are magnetically 
dominated, turbulence can be suppressed and there may 
be sufficient time for radiative levitation to operate. Fu- 
ture radiation-MHD simulations of accretion disks with 
AGN-like parameters are required to assess whether the 
surface layers are likely to be well mixed. 

6. CONCLUSIONS 

Given its astrophysical importance, the observational 
determination of SMBH spin is clearly of great interest. 
The study of X-ray reflection features from the innermost 
regions of SMBH accretion disks is, at the current time, 
the most promising technique for measuring SMBH spin 
and is a prime focus of future X-ray observatories such 
as Astro-H and LOFT. Thus, now is the time to test 
and refine our methodologies for extracting robust and 
reliable black hole spins from the X-ray data. 

The discrepant published values for the SMBH spin 
in NGC 3783 provide an ideal opportunity to examine 
and improve our methodology Here, we identify sev- 
eral factors that have led to this discrepancy. For AGN 
that are strongly affected by a warm absorber, computa- 
tional practicalities dictate that we describe the (possibly 
multiple) warm absorber zone(s) using tabulated grids 
of pre-calculated photoionization models. We point out 
the importance of using grids that adequately sample 
parameter space, especially in the ionization parameter 
£ dimension. We propose a 5-sample-per-decade crite- 
rion and show that the linear interpolation of warm ab- 
sorber grids have significant inaccuracies if sampled at a 
significantly lower density. Most importantly, however, 
we identify a statistical coupling between the iron abun- 
dance of the accretion disk and the inferred black hole 
spin. We show that one must allow for the possibility 
of non-solar abundances. Putting everything together, 
we find that the data for NGC3783 require a moderately 
supersolar iron abundance (Z ~ 2 — 4Z ) and a rapidly 
rotating black hole (a > 0.89). 

While the abundance enhancement required for 
NGC 3783 is modest and within the realm of metallic- 
ity estimates from BLR studies, we note that the X-ray 
analysis of some other AGN have found extreme iron 
abundances. This prompts us to consider mechanisms 



that can enhance the observed abundances above and 
beyond the true macroscale abundances of the circum- 
nuclear region. We give a plausibility argument to show 
that radiative levitation, the radiation-driven diffusion of 
iron towards the disk photosphere, may operate in the in- 
nermost regions of high accretion rate AGN disks. More 
theoretical work is required to assess whether radiative 
transfer effects and/or turbulent mixing of the photo- 
spheric regions of the disk will render radiative levita- 
tion ineffective. However, even the current calculation 
makes predictions that can guide future observational 
study. The clearest prediction is that there should be 
a positive correlation between the Eddington ratio of an 
AGN and the iron abundance inferred from X-ray reflec- 
tion spectroscopy of the innermost disk. This may show 
up in detailed studies of disk reflection across a sample 
of AGN, or detailed studies of a single source that is un- 
dergoing accretion rate changes. Indeed, if we find an 
individual source that shows clear and compelling evi- 
dence for a change in its inner disk iron abundance, we 
will be forced into accepting that some mechanism (such 
as radiative levitation) is modulating the observed iron 
abundance away from that of the general environment. 

It is also interesting to consider the wider implications 
for extreme iron enhancements due to radiative levita- 
tion. Since this mechanism only operates in the very 
central regions of the accretion disk (tdriftAdyn oc r 3 / 2 ), 
it is likely to be ineffective at modifying the abundance 
of accretion disk winds that originate from larger radii 
(which manifest themselves as X-ray or UV absorbers). 
However, the surface regions of the innermost disk are 
very likely the source of the "funnel flow" which bounds 
the base regions of Blandford-Znajek powered relativis- 
tic jets (Blandford & Znajek 1977; McKinney & Narayan 
2007). Thus, especially in their base regions, Blandford- 
Znajek jets may entrain plasma that has been highly en- 
riched in iron with potentially important consequences 
for high-energy hadronic processes within jets as well 
as the composition of high-energy cosmic rays that may 
originate from AGN jets. 
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